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ABSTRACT 

This study contains the analysis of the costs associ- 
ated with the printing operation of the messages received 
from the Naval Communications Fleet Multi-Channel Broad- 
cast System by a ship equipped with the Multiplex General 
Address Reading Device (GARD) . A computer simulation is 
used to compare the average message delays for three dif- 
ferent models of printers in varying quantities. Messages 
are accounted for by all four precedences. Costs of mes- 
sage delay and initial purchase price are used to compare 
the alternatives. The results of the study show that two 
high speed printers such as the TT— 624 is the optimal con- 
figuration . 



2 



TABLE OF CONTENTS 



I. INTRODUCTION AND STATEMENT OF THE PROBLEM. . . 6 

II. APPROACH TO THE PROBLEM 8 

A. THE QUEUEING MODEL 8 

B. BREAKDOWNS 11 

C. COST MODEL . . . 16 

III. SOURCE OF DATA AND ANALYSIS 18 

A. MESSAGE INTERARRIVAL TIME DISTRIBUTION 

AND RATE 18 

B. MESSAGE LENGTH DISTRIBUTION 20 

C. COST INPUTS 21 

D. MISCELLANEOUS INPUTS 2 3 

IV. THE MODEL 2 5 

A. VERIFICATION OF THE MODEL 2 5 

B. RESULTS 2 7 

C. DISCUSSION OF THE RESULTS ' . . . 27 

V. CONCLUSIONS. 31 

APPENDIX A 33 

COMPUTER PROGRAM 38 

LIST OF REFERENCES 4 2 

INITIAL DISTRIBUTION LIST 43 

DD FORM 1473 45 



3 



LIST OF TABLES 



I. Miscellaneous Inputs 24 

II. Summary of Computer Runs to Include 

Breakdown Runs 2 8 

III. Tabulated Final Results 29 



4 



LIST OF FIGURES 



1. Typical Breakdown 13 

2. Cumulative Distribution Function 

of Message Lengths 22 



5 



INTRODUCTION AND STATEMENT OF THE PROBLEM 



I . 

This study concerns the costs associated with the 
printing operation of the messages received from the 
Naval Communications Fleet Multi-Channel Broadcast System 
by a ship equipped with the Multiplex General Address 
Reading Device (GARD) . The Multiplex CARD unit performs 
the function of gathering all messages that are received 
by the various channels and then feeding the messages to 
a printer or printers. The GARD unit also performs the 
function of "weeding out" messages that are not addressed 
to that particular ship by reading the coded heading and 
then only sending those messages to the printer that are 
addressed to that ship. This aspect of the GARD was its 
primary reason for development. However, the GARD also 
eliminates the need for exactly as many printers as chan- 
nels that the ship is monitoring. In the case of a des- 
troyer, this is four channels and four printers. The 
problem, then, becomes a queueing problem with the GARD 
acting as the waiting line area and the printer or print- 
ers acting as the service channels in the model. The ob- 
jective is to minimize cost, where cost includes a charge 
for message hours lost, over a range of the number of 
printers and also three types of printers that are cur- 
rently in stock or are being evaluated for future purchase. 
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The three types of printers that were investigated in 
this study were (1) the Model 28— the present teletype 
used in most ships with a speed of 100 words per minute 
(wpm) , (2) the Diablo Model 1200 Hytype I Serial printer 

with a speed of 360 wpm, and (3) the TT— 624 produced by 
Data Products with a speed of 3000 wpm. The first was 
looked at for obvious reasons. The Diablo printer has 
been tested to some degree by the Naval Electronics 
Laboratory Center (NELC) , in San Diego, California, and 
is considered a candidate for future purchase by the 
Navy. The TT— 624 has undergone quite extensive testing 
as a high speed printer in the laboratory as well as at 
sea . 

The problem is to measure in dollars the various 
costs associated with the different printers and differ- 
ent numbers of printers. This study is a model to cal- 
culate these costs and to make a recommendation on an 
optimal type and number of printers. The costs to be 
compared are those incurred by the system to include only 
the printing operation and the cost of the queue that the 
GARD system builds up. 
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II. 



APPROACH TO THE PROBLEM 



In order to determine the optimal type and number of 
printers, the costs associated with the various setups 
must first be determined. Since the message traffic flow 
through the system behaves like a waiting line, the appro- 
priate type of model which best fits the situation is a 
queueing model. 

A. THE QUEUEING MODEL 

To measure the costs connected with a queue the 
average number of customers in that queue must be measured. 
In this case, the customers are messages. It seems 
reasonable that the relative cost of a customer in a queue 
is proportional to the importance of the customer. This 
importance of the customer is measured in communications 
by assigning a precedence or priority to the message. A 
flash message is the most important, then immediate, 
priority, and finally routine. The model must somehow 
measure the average number of messages of each precedence 
in the queue. In queueing theory this type of queue is 
called a priority queue. Since the alternatives under 
investigation include up to a maximum of four printers, 
the theoretical problem is thrown into the category of a 
multi-channel, or s server queueing model. Also, since 
the arrival rates and service rates for the messages of 
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different precedence are not necessarily equal, the 
model selected from the theoretical library of queueing 
models must maintain provisions for these characteris- 
tics . To summarize, the model to be used must be a 
multi-channel priority queue with provisions for unequal 
arrival rates and unequal service rates among the differ- 
ent precedences. Because this particular combination of 
characteristics makes the derivation extremely complica- 
ted, the solution to this problem is not found in the 
literature. For example, Cobham [Ref. 2] derived results 
for a single server priority queue with Poisson arrivals 
and a general but unequal service time distribution for 
the different priorities. He also gave results for a 
multi— server priority system with Poisson arrivals and 
the service times exponentially distributed with the ser- 
vice rate equal for all priorities. 

The author did some work on a simpler model that in- 
cluded two priorities, one channel, Poisson arrivals, and 
exponential service times. These results, found in 
Appendix A, are not new but the approach is different. 
Attempts to solve the same problem with two servers could 
not be solved using the same technique. It was obvious 
that this type of problem becomes extremely cumbersome 
and uniquely complicated whenever the number of channels 
exceeds one. 
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Since none of the theoretical models fit the problem, 
a computer simulation programmed in the General Purpose 
Simulation System/360 (GPSS) language was written and is 
included. The basic idea in the simulation is to use 
what the language calls a storage location to simulate 
the printing operation. The user can then merely set the 
storage capacity to any integer to correspond to the num- 
ber of printers being simulated. A transaction (message) 
trying to enter the storage location is permitted to do 
so only if the storage location is not full (not all 
printers are printing) . Each transaction occupies one 
unit of storage (utilizes one printer) . 

Messages are generated by the program according to 
the distribution specified by the user. Each precedence 
is generated separately but according to a certain 
theoretical mix. A message is then assigned a message 
length according to a distribution specified by the user. 
Flash messages receive a special treatment at this point. 
A 25 word flash message heading is added to the message 
length of a flash message to simulate an advance warning 
that is placed at the head of all flash messages. Next 
the message length of messages of all precedences is con- 
verted to a time that the message will require to be 
printed. Then the message joins a queue that is waiting 
for access. to the printer or printers. The above steps 
are performed for each precedence of message. It should 
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be noted that all the steps before the message enters the 
queue require no clock time in the simulation but the 
list of events is merely given to show the reader the se- 
quence of these events. Queue statistics for each prece- 
dence queue are maintained automatically by the GPSS 
program. At this point all messages from the four queues 
vie for the storage location (printer or printers). Mes- 
sages are then printed by precedence. If two messages of 
different precedence are competing for the printer at the 
same time, the one with the higher precedence will gain 
access when a printer becomes available. The queue disci- 
pline is not preemptive, that is to say, a message does 
not knock a lower priority message from the storage loca- 
tion (printer or printers). Immediately upon completion 
of printing, the message leaves the queue in which it was 
originally a member. This means that the queue statistics 
that are maintained actually include the time spent wait- 
ing to print, if any, and the printing time. The desired 
average number of messages in the printing queue or being 
printed is output as average queue contents for the queue 
of each precedence. 

B. BREAKDOWNS 

A convenient way to simulate a breakdown or a failure 
of a printer would be to merely introduce into the queue- 
ing model a transaction (in this case a breakdown instead 
of a message) w T ith a precedence higher than any other 
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customer. The user could specify the inter-arrival time 
distribution and the mean time between arrivals would be 
the Mean Time Between Failures (MTBF) of the printer. 

The service time distribution would be the repair time 
distribution with its mean time as the Mean Time to Repair 
(MTTR) of the printer. However, because the MTBF of the 
printers investigated is so long in comparison to other 
portions of the simulation, this procedure was not used. 

The simulation uses a time step routine and to get accu- 
rate queue statistics, some of the printers require the 
time step to be tenths of seconds. To simulate an adequate 
number of breakdowns would be a very inefficient use of 
the computer facility. 

Because the ratio of the MTBF to MTTR is at least 400 
to one for all alternatives, the following approach was 
taken for the alternatives involving only one printer. 

The number of messages in a queue before, during, and 
after a breakdown plotted versus time might look something 
like Figure 1. Then let: 

x = stochastically determined time to repair printer. 

t = time to restore queue to equilibrium and set to 

t = kx where k is some positive real number to in- 
sure that the jam caused by the breakdown has 
cleared . 

= average number of messages of precedence j in the 
system in equilibrium with one printer and no breakdowns. 
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Figure 1. Typical Breakdown 



N q = average number of messages of precedence j in the 
system during the time t. This is not an equili- 
brium condition since equilibrium does not exist 
during a breakdown. 

X = breakdown rate or 1/MTBF 
13 

Now for one breakdown, the "wasted message time" caused by 

the breakdown is N„.t — N. .t. Because the X^ is so small 

Oj lg B 

the following expression for the average number of messages 
of precedence j in the system for one printer including an 
average breakdown is : 

N, • = N, . + X_ (N rt . t - N, . t) 

1 ] I 3 B Oj 1 ] 

A breakdown package in the program when inserted simulates 
breakdowns in succession. As determined by the user, the 
coefficient k determines the time to insure that the jam 
of messages caused by the breakdown has cleared. At this 
time another breakdown is simulated. Then if n breakdowns 
are simulated in succession and T, obtained from the com- 
puter output, is the total time of the simulation run, set 
t = T/n and the equation above becomes 

Nij = Nlj + >. B (N oj T/n - M ljT /n) = N 1;j + V /n(N 0j " N lj> 

When i ( i > 1 ) printers are being simulated, such a dis- 
aster does not occur when one printer fails. The queues 
will merely tend to increase during the breakdown to a 
higher equilibrium level. It should be pointed out that 
for those arrival rates of messages considered in this 
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study, all three models of printers were capable of han- 
dling the load with only one printer. This is to say 
that an equilibrium exists for all models even with just 
one printer when there are no breakdowns. 

Now if: 



N^_. = average number of messages of precedence j in the 
system with i printers working in equilibrium and 
with no breakdowns. 

A = repair rate or 1/MTTR 

Then A / (A + A ) is the expected fraction of time that a 

D D K 

printer will be broken down. Let p = A r) /(A_ 1 + A^) and if 

D D K 




ever there are one or less printers operating and P^ = 
i k 

(jJ p (1— p) " is the long run probability that there are 

k printers out of i working, then 



i 




k=2 



P, N. . 
k kj 



+ (1 



i 



E V °ij 

k=2 



and is the average number of messages of priority j 

in the system when there are i ( i > 1 ) printers being simu- 
lated. Although is not readily computable it is 

evident that: 






< N 



lj 



This is true because N, . was calculated 

ID 

a breakdown of only one printer. This 

repair rate during the breakdown is A_, 

K 



while simulating 
implies that the 
or the expected 
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repair time is 1/A . When looking at the portion of time 

that all i printers are failed in the calculation of , 

the repair rate is iX R or the expected repair time is only 

1/iA^. This means that Q, . would be less than N, . because 
R l] lj 

a breakdown would be repaired more quickly on the average. 

These heuristic results complement the computer simu- 
lation to save computer time that would have been mostly 
nonproductive waiting for another breakdown to occur. 



C. COST MODEL 

One of the most difficult tasks in any model involving 
costs is determining which costs are relevant and what 
values to put on the costs that are considered relevant. 

In this model surely the costs involved with a delayed mes- 
sage are to be included. The initial cost of the printers 
is a relevant cost. Maintenance costs are another item 
that should be investigated. However, when the Navy pur- 
chases an item such as a printer in a large quantity, a 
maintenance agreement is usually included in the purchase 
price to include repair parts. Since maintenance personnel 
available for repair work would probably be the same for 
all alternatives, these costs were omitted. For the pur- 
poses of this study, then, maintenance costs were not used 
in the computation of the total cost. 

If CL is the cost in dollars per hour of delaying a 
message of precedence j and N^j is the average number of 
messages of precedence j in the system as derived in 
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Section II. B., then total cost, TC, of a selected alterna- 
tive using i printers is given by: 



4 



TC . = > v N. .C . + iP/DL 
i^i]] 

j=l 

where P is the price in dollars of the printer involved 

and DL is the design life of the printer in hours. If the 

concept of present value and discounting of costs are used 

the formula becomes: 

4 



- E 

j-i 



N. .C . 
11 1 



+ iP- 



m 



1— e 



~(m) (DL) 



where m is the annual discount rate and set at 0.10 by the 
DOD. The latter form of the cost equation considering pres- 
ent value was used in the model. 
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III. SOURCE OF DATA AND ANALYSIS 



There were basically three inputs to the model that 
were required: (1) message interarrival time distribu- 

tion and rate of arrivals, (2) the message length distri- 
bution, and (3) cost or value of reducing the writer— to— 
reader time of a Naval message. 

A. MESSAGE INTERARRIVAL TIME DISTRIBUTION AND RATE 

It would be natural to suspect that messages arriv- 
ing at a NAVCOMMSTA do so according to a Poisson process. 
Oelmann [Ref. 4] reported that this was in fact true for 
data taken during the period 5—8 March 1971 at NAVCOMMSTA 
San Francisco. lie used the Chi— square goodness— of— fit 
test to verify this result. Then it is assumed that the 
fraction of this traffic that is transmitted to a partic- 
ular ship is selected randomly from all incoming traffic 
to the NAVCOMMSTA. This is to say, each message that 
comes into the NAVCOMMSTA has a common probability of be- 
ing destined for a particular ship. This assumption 
guarantees that the traffic transmitted to a particular 
ship is still a Poisson process. The probability that 
any message is destined for a particular ship must be in- 
dependent of the destination of any other message sent 
before or after. Since there was no evidence to the con- 
trary, this assumption that messages destined for a 
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particular ship arrive at the transmitter according to a 
Poisson process was accepted. 

Given that the messages for a certain ship arrive at 
the transmitter according to a Poisson process does not 
guarantee that the output of the transmitter is a Poisson 
process, however . This would be the case if service times 
are distributed according to the exponential distribution 
(see Burke [Ref. 1]). In this case the service time (time 
to transmit) is not exponential but is distributed accord- 
ing to the message length (see Section III.B.). It should 
be clear, however, that if the traffic is light, the out- 
put should approach a Poisson process regardless of the 
transmission time distribution. The fact that traffic 
from four channels operating more or less independently is 
lumped at the GARD unit also assists in making the inter- 
arrival times of messages at the GARD unit look like an 
exponential distribution. 

A simplified transmission system consisting of four 
independent channels was simulated as described above and 
the interarrival times at the receiver were tabulated and 
a Chi-square goodness— of— fit test was accomplished. These 
results did show that the distribution of interarrival 
times at the receiver and the GARD unit to be from an ex- 
ponential distribution. 

The rate of message arrivals was determined from data 
taken aboard the USS Jouett, a DDG, in October, 1969 [Ref. 
5] . Although traffic intensity varied somewhat during the 
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day, an average rate of arrival was found to be 1/4 mes- 
sage per minute. Another source of data was that taken 
by Oelmann [Ref. 4] at NAVCOMMSTA San Francisco in March 
1971. The rate indicated there was approximately 1/2 
message per minute to all classes of ships. Since this 
rate also varied considerably from day to day and during 
the day, it v/as felt that a precise rate would be hard to 
determine. Therefore the rate of 1/4 message per minute 
was selected and thought of as a fairly heavy rate for 
this class of ship. It should be noted that although 
this may be a fairly heavy rate for an average value, it 
is probably quite low for some heavy traffic periods that 
all ships encounter. 

B. MESSAGE LENGTH DISTRIBUTION 

The message length distribution was required because 
the length of a message is what determines that time that 
a message will take to be transmitted or printed. As 
stated earlier, in the queueing model this turns out to be 
the service time of the transmitter or printer. It is 
clear that the service time is scaled by the speed of the 
printer being utilized but the message length alone deter- 
mines the service time of any model or type of printer. 

The speed of the transmitter is held constant at 100 words 
per minute, the speed at which the Fleet Broadcast is 
actually transmitted. 

The message length distribution was found by Oelmann 
[Ref. 4] to be a right shifted hyperexponential 
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distribution. The source of the data was a sample of mes- 
sages taken at NAVCOMMSTA San Francisco on 5 October 1971. 
Another source of data was the Naval Electronics Labora- 
tory Center in its World Wide Traffic Demand Study data 
bank. Those data were taken on 8 October 1970 in various 
parts of the world. The determination of the distribution 
as a right shifted hyperexponential was confirmed by the 
Chi— square and Kolmogorov— Smirnov (K— S) goodness— of— fit 
tests. The reason that the distribution is right shifted 
can be explained by the message heading that appears on 
all messages regardless of precedence. This heading con- 
tains a data— time group, an originator, and the addressee. 
Figure 2 shows the cumulative distribution function (CDF) 
for the four precedences of messages. 

Both sets of data were accepted by the goodness— of— 
fit as originating from a right shifted hyperexponential 
but the data taken at NAVCOMMSTA San Francisco was used 
to determine the parameters for the distribution in the 
simulation . 

C. COST INPUTS 

Another input to the basic cost model was the price 
of time that a message was delayed. In the model it was 
the price of time that a message spent in the buffer be- 
fore being printed and the time being printed. Washburn 
[Ref. 6] performed two experiments to determine the price 
of time of a Navy message while working at the Naval 
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Figure 2. Cumulative Distribution Function of Message Lengths (in words). 



Communications Command in December 1972. The most rele- 
vant experiment involved using the Delphi technique to 
survey a group of senior Navy Communications personnel. 

The technique involves feedback to the person completing 
the questionnaire in between repeated personal estimates. 
After the first survey, all participants are shown the 
results of that survey before being surveyed again. Thus, 
a member's estimate will be based upon his knowledge of 
all other participant's first response. The quantity of 
interest in the experiment was the price of time for 
Naval messages. This can be defined as the amount of 
money the Navy should be willing to pay to shorten the 
writer— to— reader time of a message by one hour. All esti- 
mates were instructed to be of average value. The results 
of the third and final survey taken of ten senior communi- 
cations personnel gave median values of $1 per hour for a 
routine message, $10 per hour for a priority message, $100 
per hour for an immediate message, and $1000 per hour for 
a flash message. 

It should be noted that although the estimates varied 
widely, the median values of price had reached the final 
values after only two surveys. The third survey merely 
lowered the variance within the sample. 

D. MISCELLANEOUS INPUTS 

Mean Time Between Failures (MTBF) and Mean Time To Re- 
pair (MTTR) estimates for the three models of printers 
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were either obtained from manufacturer's data or from ex- 



perienced technicians at the Naval Electronics Laboratory 
Center at San Diego. Table I summarizes the values ob- 
tained and used in the model. The costs of the printers 
were also obtained from NELC and are included in Table I . 
Department of the Navy has set the design life of all 
minor communications equipment such as printers at eight 
years . 

The fraction of messages by precedence that were used 
in the model were 1% flash messages, 15% immediate mes- 
sages, 54% priority messages, and 30% routine messages. 
These figures were taken from actual data taken aboard the 
USS Jouett in October 1969 and substantiated by other 
studies to be meaningful figures (see [Ref. 5]). 









Characteristic 




Model 


MTBF 


MTTR 


Purchase 

Price 


- Design 
Life 


Model 28 


800 hrs 


2 


hrs 


$ 1500 


8 yrs 


Diablo 


1500 hrs 


3 


hrs 


$ 1900 


8 yrs 


TT— 624 


3000 hrs 


3 


hrs 


$15000 


8 yrs 



TABLE I. 

Miscellaneous Inputs 
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IV. THE MODEL 



The GPSS program as outlined in Section II was uti- 
lized to provide the average queue contents statistics 
that were required for the cost model. - When the break- 
down package was used, field A of the TERMINATE block 
for all the messages was set to zero so that a specific 
number of breakdowns would terminate the simulation run 
rather than a combination of breakdowns and messages. 
Field A of the TERMINATE block is the amount by which 
field A of the START block is decremented for each trans- 
action that comes into the TERMINATE block. The simula- 
tion terminates when field A of the START block is less 
than or equal to zero. 



A. VERIFICATION OF THE MODEL 

Cobham [Ref. 2] gives the following results for an 
M/G/l queue with r priorities of customers. Let S^ be 
the service time of a customer of precedence j . Let 
Fj (t) be the CDF of the service time distribution for 
customers of the jth precedence and A ^ is the arrival 
rate for the customers of the jth precedence. Let r be 
the highest priority and 1 be the lowest precedence. 
Then 

r r 

A = £ A. and F ( t ) = £ (A ./A) F . (t) . 



j = l 



j = l 
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Then 



dF(t)/dt = X (A./A)f. (t) 
j=l 3 3 



where f ^ (t) is the probability density function of 



t 2 dF(t) = X (A ./A) t 2 f . (t)dt 
j = l 3 3 



/ t 2 dF(t) = / X (A ./A) t 2 f . (t) dt 
> o j=l 3 3 



= (1/A) X A • / t 2 f . (t)dt 

j=l 3 o 3 



but / t 2 f.(t) = Var (S . ) + (E(S .)) 2 
o 3 11 

Then is the average waiting time in the queue of a mes- 
sage of precedence j and is given by the expression: 



W. = 
1 



(1/2) A / t dF ( t ) 
o 



j A 



i-E 



r+l—k 



k=l 



r+l—k 



31 X r+l-k 

1- X y 



k=l 



r+l—k 



where 1/y^ is the expected value of the service time or 
E ( S j ) of a customer of precedence j. Then the expected 
time in the system of a customer with precedence j is 
given by 

WS . = W . + E (S . ) 

1 1 1 



A comparison using these results and 'the results from 
the simulation for the case of one Model 28 printer was 
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performed. The results of the comparison showed that the 
expected waiting time in the system from the simulation 
was within 5% of the theoretical results for the three 
highest precedences and only slightly more for the routine 
messages . 

B. RESULTS 

All results from the simulation runs are. shown in 

Table II. A 95% confidence that all values except the 

breakdown statistics are within 5% of their true value 

was attained. For the N_ . statistics a 90% confidence 

03 

that the results were within 10% of their true value was 
attained. Table III shows the average number of messages 
in the system to include the provisions for breakdown of 
the printers. The table also shows final cost figures 
for message delay, cost of printers and the total cost 
for each alternative. 

C. DISCUSSION OF THE RESULTS 

It is evident from the results that the optimum num- 
ber of printers is two and that the TT— 624 is the most 
economical model. Even if the costs of a delayed message 
are cut by a factor of ten, the same alternative is only 
slightly less economical than two Diablo printers. The 
primary driving parameter in the model is the speed of 
the printers. The slower printers waste too much time 
merely printing even when there are adequate numbers of 
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TABLE II. Summary of Computer Runs to Include Breakdown Runs 
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TABLE III. Tabulated Final Results 



them to virtually eliminate a queue in the buffer. At 
times during the day when traffic is heavier, the choice 
of the TT— 624 printer will be even a better one than dur- 
ing lighter traffic. Some runs were made at a rate of one 
message per minute and these results substantiate the 
above statement. It should be noted, however, that the 
alternatives of one or two Model 28 printers could not be 
simulated at this high rate because it is theoretically 
impossible for this model to handle a one message per 
minute arrival rate with only one or two servers. 
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V. CONCLUSIONS 



The conclusion that can be drawn from this study is 
that two TT— 624 printers is the most economic choice 
from among those alternatives investigated. It should 
be noted at this point that a lot of the stated assump- 
tions in the model are unimportant or irrelevant in the 
final analysis. The model points out a rather common 
sense result — one high speed printer is sufficient but 
that a spare is required when the primary printer breaks 
down . 

An alternative that was not investigated because of 
programming difficulties is a mixture of two different 
models of printers. An example is a high speed printer 
to act as the primary printer and handle traffic during 
the majority of the time and a Model 28, which is inex- 
pensive, to maintain a minimum queue during the period 
when the high speed printer is broken down. This would 
involve a slightly more complex GARD unit and would in- 
volve extra repair parts so it is not immediately 
apparent if the advantages would outweigh the disadvan- 
tages . 

Space on a ship such as a destroyer is also an im- 
portant factor and has a cost associated with it. The 
most economic alternative involves the saving of some 
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space as compared to the four Model 28 teletypes that 
are presently installed in a destroyer. It was not con- 
sidered in the study but it is pointed out here to show 
another advantage of the alternative chosen. 

It might be argued that the Model 28 teletypes that 
are currently in service are free but even if this is 
the case, the Model 28 teletype is still inferior because 
of its slow speed. It should be pointed out, however, 
that costs such as retraining costs for a high speed 
printer maintenance course are not included in the study. 
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APPENDIX A 



SOME THEORETICAL RESULTS 



In the initial stages of this study, the author at- 
tempted to derive the theoretical results which were not 
found in the literature. A simple problem for which the 
results were known was first used to test a solution 
technique. Given an M/M/1 queue with two priorities of 
customers, customers of the higher priority have inter- 
arrival times distributed EXP(X 2 ) and service times dis- 
tributed EXP (y 2 ) ■ Lower priority customers have 
interarrival times distributed EXP(X^) and service times 
distributed EXP(y^). If the queue discipline is first— in- 
first— out within each priority class, the higher priority 
customers preempt the ones of lower priority, and P(i,j) 
is the steady state probability that there are i high 
priority and j low priority customers in the system, then 
the balance equations become: 



(X 1 +X 2 )P(0,0) =y 2 P(l,0) +y 1 P(0,l) 

(X 1 +X 2 +y 1 )P (0 , j ) =A 1 P (0 , j-1) +y 1 P (0 , j+1) +y 2 P (1 , j ) for j >1 
(X 1 +X 2 +y 2 ) P (i , j ) =X- L P (i , j — 1 ) +X 2 P ( i— 1 , j ) +P 2 P ( i+1 ' 3 ) for i-1 1 > 



j>0 and P(i,— 1)=0 and 



E E 

i=0 j =0 



P (i , j ) =1 



Then let Q (i , z„ ) = P(i,j)z^ and R(z ,z„)=^ Q(i,z„)zj 

j=0 i=0 

where CKz^,z 2 £l 
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Then the balance equations can be worked into the form 



M 1 Z 1 ( Z 2 _1 ) P (°' °) + [y 2 Z 2 (z 2 ~l) ]Q (0 , z 2 ) 

_X 1 Z 1 Z 2 + Z 2 (A 2 Z 1 + X 1 Z 1 +1J 2 Z 1 _A 2 Z 1 _M 2 ) 



At this point it would be well to point out that although 
this expression is still in terms of the two unknowns 
P(0,0) and Q(0,z 9 ), the fact that 

Z co 00 

X) ^ P (i, j )=1 

i=0 j=G 



(or equivalently that R(l,l) = l)has not been used. After 
this transformation to the (z^,z 2 ) space it is now hoped 
that by approaching the point (1,1) from different direc- 
tions that more information about the above equation can 
be found to solve for these unknowns. 



Then if z 2 = f(z^) = 1+a (z^~ 1) +b (z^— 1) and f (1) = 1 



f 1 (z^) = a+2b(z^— 1) 
f " (z JL ) = 2b 



f ’ (1) = a 

f " (1) = 2b 
0 



Now let R(z.,z 2 ) = R(z^,f(z^)) and lim R(z.,f(z^)) = ^ 

z^->l 

so using L' Hospital's rule yields 



R(l, f (1) ) = 



(m 1 P(0,0)-u 1 Q(0,1) ]a+M 2 Q(0,l) 



t — a+y 2 -X 2 



where Q(0,1) = Q(0,f(z 1 )) 



z l=! 



but R ( 1 , f ( 1 ) ) =1 for all "a". Therefore it is deduced that 
Q (0 , 1) = 1-A 2 /m 2 = l-p 2 and then P(0,0) = 1-X ± / \i ± -X 2 / 2 = 

1— p^— P 2 . Then if a = (y 2 — X 2 )/X^, L'Hospital's rule can be 

used again. Again equating terms in the numerator and 
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denominator that do not have b as a coefficient because 
R(l,f(l))-1 for all b, an expression is obtained for 



Q' (0,f (z 1 ) ) 



Z l =1 



It is Q' ( 0 , f ( 1 ) ) = 



p l [ (1 "' P 2 )2+U 2 A 3 )/,J 2 ] 

(1-P 2 ) (l-p 1 -p 2 ) 



Now 



R(z lf l) = 



p 2 (z 1 1 ) Q ( 0 , 1 ) _ Q ( 0 , 1 ) 



A 2 Z 1 +M 2 Z 1 A 2 Z 1 U 2 1 P 2 Z 1 



Then the expected number of high priority customers in 



the system is d/dz^R(z^,l) 



Z l =1 



p 2 A 2 



which is a well known result from queueing theory consider- 
ing that the high priority customers act as if they are 
alone in the system and are not influenced by the low 
priority customers. 



Now 



R(1,z 2 ) = 



p 2 ? ( 0 , 0 ) — y 2 Q ( 0 , z 2 ) 
“ X 2 Z 2 



Then the expected number of low priority customers in the 
system is 



d¥ 2 R(1 ' z 2 ) 



„ = y— (Q ' ( 0 , 1 ) +P ( 0 , 0 ) — Q (0,1)) 

z 2 -l A x 

P 1 P 1 
- y^(Q ' (0 , l)-p, ) = ^'(0,1) -1 

A 1 1 A 1 



and inserting the expression for Q' (0,1) that was derived 
above gives the result 



dz 2 R(1,z 2 ) 



= A 1 P 2 /P 2 +P 1~ P 1 P 2 
Z 2 =1 (l-p 2 ) (l- Pl -p 2 ) 



This agrees with results found by Jaiswal [Ref. 3]. 
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Using the same technique as in the above problem to 
solve an M/M/2 queue with the same two classes of customers 



as in the above problem led to a more complex expression 
for R(z^,z 2 ) that involves Q(0,z 2 ), P(0,0), and P(0,1). 

(A) Q ( 0 , z 2 ) + (B ) P (0,0) 4- (C)P(0,1) 

FU Z 1' Z 2^ - ; ' l i + ^2 +2 ’ J 2~ 2 ’ J 2 //z l~ ; ' l l Z 2 — ^2 Z 1 



where 



A= (z — 24— 
1 U 



Ul z l 



2 Z 2 



y l Z l 2u l 

U 1 +X 2 +2u 1 X 1 z 2 — ) 



+ (X 1 +X 2 +2u 2 - 



2y . 



- \ x z 2 ) 



B= (z — 24- 



^1 Z 1 

U 2 z 2 



M 1 z . 2p 

-i-i) (l-z 2 ) (X 1 + X 2+ -+) 



C ~ ( Z 2~ 1 } { —kr + P 2 z l “ 2i1 2 “ lJ l Z l } 



Again setting z 2 = f(z^) = l+a(z^— 1) 4- b(z^— 1) and using 

L' Hospital's rule gives a relation between P(0,0) and 

P(1,0) as well as the value of Q(0,1). Again selecting 

"a" properly allows a second application of L' Hospital's 

rule but this time no additional information is found about 

the relation between P(0,0) and P(1,0) when the terms with 

a coefficient of b in the numerator and the denominator 

are equated. Again Q' (0,1) appears but its solution is in 

terms of P(0,0). Then to take — R(l,z„) yields an expres— 

dz 2 Z 

sion containing Q 1 (0,1) and P(0,0). Therefore the final 
expression for the expected number of low priority custom- 
ers in the system is still in terms of the unknown quantity 
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P(0,0) and cannot be solved. 

E (Number of low priority customers in the system) 



~ dz 2 R(1,Z 2 


(A) P (0, 0) + B 
z 2 =i Pi 


where 





A = (2+p 2 ) (2y 2 +2y 1 +A 1 +X 2 ) / (l-y 1 a/y 2 ~2y 1 a-2y 2 ) 

B = 2y 1 /y 2 -p 2 -A 1 /y 2 -2+2A 1 -4y 2 /y 1 -2ap 1 + [2-p 2 ] [p 2 y x (2-a) 
+y 1 a (y 1 /y 2 -2-A 1 /y 2 ) +A X J -(2y 1 ) (2-p 1 ~p 2 ) (2+p 2 )/ (l-y 1 a/y 2 



— 2y^a+2y 2 ) 




where a = (2y 2 ‘ 


-x 2 )/q 



One possible use of an equation such as this is that 
it is possible that the quantity P(0,0) may be easier to 
measure in a particular case than the number of low prior- 
ity customers in the system. Here this equation could be 
used to find the expected value that was desired. 
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